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ABSTRACT 


We  report  concurrent  radio  and  gamma-ray  observations  of  PSR  B1509  —  58  carried  out 
by  the  Parkes  Radio  Telescope  and  by  the  Burst  And  Transient  Source  Experiment 
(BATSE)  and  the  Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  the 
Compton  Gamma  Ray  Observatory  ( CGRO ).  Gamma-ray  light  curves  fitted  at  several 
energies  between  ~  20-500  keV  yield  a  phase  offset  with  respect  to  the  radio  pulse 
that  is  independent  of  energy,  with  an  average  value  0.32  ±  0.02.  Although  this  value 
is  larger  by  0.07  than  that  reported  by  Kawai  et  ah,  the  difference  is  not  statistically 
significant  (only  ~  2<r)  when  account  is  taken  of  the  uncertainty  associated  with  their 
result.  We  briefly  discuss  the  possibility  that  the  energy-independence  of  the  gamma- 
ray  pulse  phase  is  a  signature  of  non-thermal  radiation  in  the  X-ray /gamma- ray  range 
and  the  suggestion  of  a  dependence  of  pulsar  radio— gamma-ray  phase  offset  on  pulse 
period. 

Subject  headings:  gamma  rays:  general  —  pulsars  —  radio  sources:  general  —  X-rays: 
sources 


1.  INTRODUCTION 

Until  recently,  there  were  only  two  radio  pulsars,  the  Crab  and  Vela,  which  had  also  been 
detected  in  hard  X-rays/gamma-rays,  and  they  differed  both  in  spectral  shape  and  in  the 
behavior  of  the  pulse  phase  as  a  function  of  energy.  The  Crab  pulsar  spectrum  can  be  fit  by 
a  power  law  in  the  ~  0.1-200  keV  energy  range  (Knight  1982,  1983,  and  references  therein; 
Ulmer  et  al.  1993).  The  double-peaked  pulse  is  aligned  in  phase  at  all  frequencies  from  the 
radio  to  the  gamma  rays.  On  the  other  hand,  the  Vela  pulsar  X-ray  spectrum  in  the  0. 1-2.4 
keV  range  is  probably  of  thermal  origin  (Ogelman,  Finley,  &  Zimmerman  1993)  and  the 
phases  of  the  Vela  pulsar  peaks  shifts  from  the  optical  to  X-rays  to  gamma-rays  (Ogelman 
et  al.  1993;  Manchester  &  Taylor  1977  (MT77);  Lyne  and  Graham-Smith  1991  (LGS91)). 

Instruments  on  the  Compton  Gamma-Ray  Observatory  (CGRO)  have  now  detected  3  ad- 
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ditional  radio  pulsars:  PSR  B1706— 44,  PSR  B1055— 52,  and  PSR  B1509— 58  (Thompson 
et  al.  1992;  Fierro  et  al.  1993;  Wilson  et  al.  1992a).  OSSE  and  BATSE  have  not  yet  de¬ 
tected  PSR  B1706— 44  and  PSR  B1055— 52  at  hard  X-ray  energies  but  we  have  seen  PSR 
B1509— 58.  Like  the  Crab,  it  has  a  non-thermal,  power  law  spectrum  from  the  soft  to  hard 
X-rays  (Seward  &  Harnden  1982;  Seward,  Harnden,  &  Eisner  1985;  Wilson  et  al.  1992a; 
Wilson  et  al.  1993;  Matz  et  al.  1993).  However,  unlike  the  Crab  pulsar,  the  radio  pulse 
phase  differs  markedly  (~  0.3)  from  the  phase  of  the  soft  X-ray  (2-11  keV)  peak  (Kawai  et 
al.  1991  (K91)). 

In  this  paper  we  report  on  pulse  phase  observations  of  PSR  B1509— 58  which  combine  data 
from  the  BATSE  (Wilson  et  al.  1992b)  and  OSSE  (Johnson  et  al.  1993)  experiments  on 
CGRO ,  covering  ~  20-500  keV  in  energy,  along  with  nearly  contemporaneous  ~  2  GHz 
radio  observations  made  at  the  Parkes  Radio  Telescope.  By  comparing  the  hard  X-ray  and 
radio  pulse  phases  we  can  verify  the  phase  shift  measured  by  Ginga  (K91)  and  determine 
whether  the  pulse  phase  varies  over  two  orders  of  magnitude  in  photon  energy.  We  describe 
the  observations  and  analysis  procedures  in  §2,  and  present  our  results  for  PSR  B1509— 58 
in  §3.  In  §4  we  briefly  discuss  our  results  in  comparison  with  other  pulsars. 

2.  OBSERVATIONS  AND  ANALYSIS 

2.1  Radio  Observations  and  Ephemeris 

The  radio  observations  of  PSR  B1509— 58  were  made  at  the  Parkes  Radio  Telescope  over 
regular  intervals  during  the  period  MJD  47931-48896,  using  two  observing  systems  centered 
at  frequencies  of  1.52  and  2.36  GHz.  Details  of  the  observing  apparatus  and  procedures  are 
given  in  K91.  To  maintain  consistency  when  comparing  results,  we  adopt  their  dispersion 
measure  (DM)  of  253  ±  1  cm-3  pc,  derived  by  comparing  the  average  pulse  profiles  at  the 
two  observed  frequencies.  Using  these  radio  data  we  derived  an  ephemeris  for  this  pulsar; 
the  parameters  are  given  in  Table  1,  where  z/,  z/,  and  i)  are  the  pulsar  frequency  and  its  first 
and  second  derivatives. 
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The  radio  phase  is  defined  by  the  UTC  geocentric  arrival  time  of  a  fiducial  peak.  By 
convention,  the  firs4  peak  after  the  beginning  of  a  day  near  the  (temporal)  center  of  the 
radio  observations  used  to  define  the  ephemeris  is  chosen;  the  reference  time  is  referred  to 
as  t0geo-  The  ephemeris  derived  gives  the  the  arrival  time  at  the  geocenter  for  an  “infinite 
frequency”  pulse;  i.e.,  a  radio  pulse  corrected  for  DM.  To  transform  from  the  geocenter  to 
the  solar  system  barycenter  (SSB),  we  employed  a  version  of  the  standard  software  package 
TEMPO  (Taylor  &  Weisberg  1989),  in  conjunction  with  the  Jet  Propulsion  Laboratory 
DE200  solar  system  ephemeris  (epoch  J2000),  kindly  provided  by  E.  M.  Standish  (private 
communication).  The  best  fit  position  of  PSR  B1509— 58  used  for  this  solution  is  shown  in 
Table  1.  The  epoch  range  for  which  this  solution  was  derived  is  MJD  47913-48928. 

2.2  BATSE  and  OSSE  Analysis 

The  BATSE  and  OSSE  instruments  are  described  in  detail  in  Fishman  et  al.  (1989)  and 
Johnson  et  al.  (1993).  Here  we  briefly  describe  the  timing  data  that  we  analyzed  for  this 
paper.  For  both  instruments,  the  data  were  accumulated  in  relatively  broad  energy  channels 
chosen  prior  to  the  observation.  We  report  here  data  from  3  OSSE  and  7  BATSE  energy 
bands,  covering  ~  20-480  keV. 

The  BATSE  data  from  the  source-facing  detectors  were  accumulated  in  64-bin  phase  his¬ 
tograms  on  board,  obtained  by  folding  the  data  over  an  integral  number  of  “mean”  pulse 
periods  for  a  programmable  duration,  normally  8.0  s.  For  27%  of  the  exposure  reported 
here,  durations  of  130  s  have  been  used.  The  hardware  performs  the  function  of  a  multi¬ 
channel  scaler,  successively  stepping  to  the  adjacent  phase  bin  every  1  / 64th  of  the  folding 
period.  Detected  counts  are  accumulated  into  the  currently  active  phase  bin,  based  on  their 
UTC  arrival  times  at  the  spacecraft.  The  period  used  for  this  on-board  processing  was  the 
mean  observed  period  in  the  reference  frame  of  the  spacecraft,  specified  to  the  nearest  fi s. 
The  period  used  includes  the  effects  of  the  Earth’s  motion  relative  to  the  source;  a  single 
value  for  the  period  is  used  throughout  the  CGRO  orbit.  The  maximum  mean  phase  drift 
occurring  over  either  integration  time  is  less  than  0.5  phase  bins.  Data  are  recorded  sep- 
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arately  for  each  of  16  energy  channels.  These  data,  along  with  the  starting  time  of  the 
accumulation  as  measured  by  a  simultaneous  latch  of  the  CGRO  spacecraft  clock,  read  out 
with  a  precision  of  1/64  ms,  are  telemetered  to  the  ground. 

In  ground  processing,  the  SSB  arrival  time  of  the  65  bin  edges  nearest  to  the  midpoint  of 
the  accumulation  are  computed,  using  software  derived  from  TEMPO  with  the  JPL  DE200 
solar  system  ephemeris  and  the  pulsar  radio  position  given  in  Table  1.  The  counts  are  then 
summed  into  64  phase  bins  fixed  in  the  SSB,  splitting  counts  linearly  into  adjacent  bins 
based  on  the  relative  overlap  with  the  SSB  bins. 

The  OSSE  pulsar  data  for  PSR  B1509— 58  were  accumulated  as  rates  in  8  ms  time  intervals; 
each  sample  is  tagged  with  its  UTC  arrival  time  at  the  spacecraft.  After  transmission  to  the 
ground,  UTC  times  were  corrected  to  the  solar  system  barycenter  using  the  same  software 
as  BATSE,  then  epoch-folded  using  the  radio  ephemeris.  During  the  OSSE  observation,  two 
detectors  were  continuously  pointed  at  the  pulsar  and  two  alternated  source  and  background 
pointings  every  two  minutes.  Pulsar  data  were  only  collected  from  a  detector  when  it  was 
pointing  at  the  source. 

We  have  taken  considerable  care  to  ensure  that  both  BATSE  and  OSSE  (in  fact,  all  CGRO 
instruments)  use  consistent  pulsar  analysis  techniques  (including  as  much  common  software 
as  possible),  and  have  verified  the  consistency  of  our  results  by  comparing  light  curves  of 
the  Crab  pulsar. 

For  both  instruments,  the  epoch-folding  procedure  computed  the  phase  (0)  of  the  pulsar 
according  to  the  formula: 


4>{t)  —  <f>(t0)  +  Vq  (f  —  to)  +  i)q  (t  —  to)2 /2  +  Vo  (f  —  fo)3/6.  (1) 

The  radio  ephemeris  for  PSR  B1509— 08,  including  the  parameters  of  this  formula,  are  given 
in  Table  1. 
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2.3  CGRO  Timing  Accuracy  and  Stability 


The  rated  stability  and  accuracy  of  the  satellite  clock  and  the  accuracy  of  the  satellite 
position  information  combine  to  give  an  expected  absolute  timing  accuracy  of  100/is.  We 
have  used  observations  of  the  Crab  pulsar  to  verify  this  value  by  measuring  the  absolute 
timing  accuracy  and  the  long-term  clock  drift.  The  stability  of  the  clock  is  important  to  our 
analysis  because  the  BATSE  observations  of  PSR  B1509— 58  cover  a  much  longer  period  of 
time  than  those  of  OSSE  (MJD  48561-48920  for  BATSE  vs.  MJD  48910-48924  for  OSSE). 

The  absolute  accuracy  of  the  CGRO  clock  was  demonstrated  by  epoch-folding  the  OSSE 
Crab  pulsar  data  using  the  radio  ephemeris  of  Arzoumanian,  Nice,  &  Taylor  (1992).  The 
primary  peaks  of  the  gamma-ray  and  radio  pulses  agreed  to  better  than  ~  0.3  ms,  the 
typical  absolute  accuracy  limit  for  radio  pulsar  observations.  We  conclude  that  absolute 
accuracy  of  CGRO  timing  data  is  at  least  this  good. 

Limits  on  clock  drift  were  determined  by  OSSE  using  light  curves  produced  from  Crab 
observations  separated  by  three  weeks;  the  data  were  accumulated  over  periods  of  one  day 
to  one  week.  Comparison  among  these  various  light  curves  showed  no  change  to  better 
than  0.01  in  phase  over  ~  6  X  108  cycles,  demonstrating  that  any  drift  in  the  CGRO  clock 
was  negligible.  BATSE  has  observed  the  Crab  pulsar  continuously  since  the  launch  of 
CGRO.  Using  the  ephemerides  of  Arzoumanian  et  al.  (1992),  with  the  assumption  that  any 
observed  change  in  arrival  phase  of  the  observed  gamma-ray  pulse  profile  is  due  to  errors 
in  the  CGRO  clock,  a  limit  of  jlois  ms  can  be  placed  on  the  drift  of  the  clock  throughout 
the  CGRO  mission  ( jloii  ms  during  the  observations  of  PSR  B1509— 58  reported  here). 

Therefore  both  the  absolute  accuracy  and  drift  uncertainties  in  the  CGRO  clock,  determined 
independently  by  the  OSSE  and  BATSE  measurements  of  the  Crab  pulsar,  are  on  the  order 
of  0.3-0. 5  ms.  Uncertainties  of  this  magnitude  are  small  (~  0.3%)  compared  to  the  150  ms 
period  of  the  PSR  B1509— 58  and  the  X-ray/radio  phase  shift,  and  are  comparable  to  the 
statistical  errors  for  the  combined  BATSE  and  OSSE  data  (see  §3). 
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3.  RESULTS 


In  Figure  1,  we  show  the  results  of  the  analysis  described  above  in  the  form  of  three  average 
light  curves:  the  radio  data,  the  20-50  keV  BATSE  data,  and  the  best  signal-to-noise  light 
curve  for  higher  energies  (50-170  keV)  from  OSSE.  Light  curves  in  all  the  CGRO  energy 
bands  will  be  reported  elsewhere  (Wilson  et  al.  1993;  Matz  et  al.  1993). 

Although  the  pulsed  emission  is  clearly  significant  in  both  gamma-ray  light  curves,  the 
statistical  fluctuations  from  bin  to  bin  are  large  enough  that  we  cannot  accurately  measure 
the  phase  of  the  peaks  directly  from  the  data.  Instead,  we  have  fitted  the  observed  light 
curves  with  a  model  pulse  profile  in  order  to  determine  the  phase  of  the  gamma-ray  peak. 
We  used  two  approaches:  The  first  was  to  fit  the  gamma-ray  light  curves  with  a  Gaussian 
peak  plus  a  constant  amplitude  offset.  Applied  to  the  light  curves  shown  in  Figure  1,  this 
method  yielded  a  relative  gamma-ray  phase  of  0.36  ±  0.01.  The  %2  for  the  fits  were  72  and 
59  for  BATSE  and  OSSE,  respectively,  for  60  degrees  of  freedom. 

Although  these  %2  values  indicate  statistically  acceptable  fits  to  the  data,  the  Gaussian 
model  does  not  reflect  the  apparent  asymmetry  of  the  light  curve.  For  this  reason  we 
applied  a  second  model,  using  the  32-phase  bin  ~  2-11  keV  Ginga  light  curve  (K91)  as  a 
template  to  fit  the  CGRO  light  curves.  The  free  parameters  for  this  fit  were  the  phase  and 
normalization  of  the  template  plus  a  constant  DC  count  level.  Model  values  for  phase  shifts 
of  a  fraction  of  a  phase  bin  were  determined  by  linearly  interpolating  between  adjacent  bins 
in  the  Ginga  light  curve. 

For  the  light  curves  in  Fig.  1,  the  best  fit  required  a  shift  in  the  previously  reported  phase 
of  the  K91  light  curve,  relative  to  the  peak  in  the  radio  light  curve,  of  0.07  ±  0.02.  Note 
that  because  we  have  used  the  K91  light  curve  as  a  template,  their  relative  phase  of  0.25  is 
treated  as  a  reference  phase  for  the  template,  with  no  associated  uncertainty.  The  implied 
relative  phase  of  the  gamma-ray  peak  from  the  OSSE  and  BATSE  data  is  then  0.32  ±  0.02. 
For  61  degrees  of  freedom,  we  obtained  a  %2  of  84  and  53  for  BATSE  and  OSSE,  respectively. 
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Thus  this  model  is  also  a  statistically  acceptable  description  of  the  data,  which  implies  that 
there  is  no  evidence  for  a  change  in  the  light  curve  shape  from  the  soft  to  hard  X-rays. 

To  measure  the  energy  dependence  of  the  pulse  phase  we  have  accumulated  light  curves 
from  BATSE  and  OSSE  over  7  and  3  broad  energy  ranges,  respectively,  covering  the  energy 
range  ~  20-500  keV.  Then  at  each  energy  the  best-fit  phase-shift  values  were  determined 
using  the  Ginga  light  curve  as  a  fitting  template.  The  results  are  displayed  in  Figure  2. 

The  weighted-average  phase  shift  determined  from  the  combined  OSSE  and  BATSE  points 
(Fig.  2)  is  0.315  ±  0.005.  The  uncertainty  on  this  result  approaches  the  accuracy  of  the 
radio  data;  it  is  purely  statistical,  however,  and  does  not  include  the  inherent  resolution  of 
the  model  light  curve  (about  0.016  in  phase)  or  uncertainties  in  the  CGRO  clock  (~  0.003) 
or  radio  peak  position  (~  0.003).  Taking  these  sources  of  error  into  account,  we  quote  our 
best  value  for  the  phase  shift  as  0.32  ±  0.02.  Separate  fits  to  subintervals  of  the  BATSE 
observations  are  consistent  with  no  time  evolution  of  this  offset. 

4.  DISCUSSION 

The  peak  of  the  gamma-ray  pulse  of  PSR  B1509— 58  is  clearly  offset  in  phase  from  the 
peak  in  the  radio.  Our  gamma-ray  light  curves  cover  ~  20-480  keV  in  energy,  compared 
with  ~  2-11  keV  for  the  K91  measurement.  When  proper  account  is  taken  of  the  revised 
uncertainties  in  the  Ginga  position,  their  phase  offset  relative  to  the  radio  is  0.25  ±  0.03 
(K91  and  Kawai,  private  communication).  The  adjustment  of  0.07  ±  0.02  that  we  have 
determined  to  this  phase  offset  is  therefore  only  marginally  significant  (~  2cr).  It  is  clear 
from  Fig.  2  that  there  is  no  significant  variation  with  energy  over  the  ~  20-500  keV  energy 
range  spanned  by  the  CGRO  data.  Evidently  then,  PSR  B1509— 58  is  like  the  Crab  pulsar 
(cf.  LGS91;  MT77)  in  this  respect.  Unlike  the  Crab  pulsar  however,  with  its  very  small 
X-ray/radio  phase  offset  (~  0.05  or  ^  0.01,  depending  on  whether  the  precursor  or  the 
main  pulses  are  used  to  define  the  radio  reference  phase  —  see,  e.g.,  MT77  and  LGS91), 
the  relative  phase  offset  observed  for  PSR  B1509— 58  is  fairly  large. 
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It  is  difficult  to  make  a  detailed  comparison  of  this  result  with  the  gamma-ray  and  radio 
pulsars  PSR  B1706— 44  and  PSR  B1055— 52  (Thompson  et  al.  1992;  Fierro  et  al.  1993) 
because  their  gamma-ray  signals  are  weak  and  their  pulse  shapes  are  not  well  determined. 
It  does  appear,  however,  that  the  gamma-ray  peaks  of  these  pulsars  are  also  not  aligned 
with  the  radio  peaks.  Since  this  is  also  the  case  for  the  Vela  pulsar  (LGS91;  MT77),  in 
this  respect,  at  least,  the  Crab  is  the  “unusual”  object  in  the  set  of  CGRO  detected  radio 
pulsars. 

We  note  (as  K91  noted)  that  in  progressing  from  the  Crab  to  PSR  B1509— 58  the  period 
and  the  radio/gamma-ray  phase  offset  seem  to  increase  together.  There  are  a  number  of 
difficulties,  however,  which  make  this  apparent  correlation  problematic.  It  is  not  obvious 
how  to  construct  a  meaningful  definition  of  phase  shift  which  can  be  uniformly  applied  to 
pulsars  with  markedly  different  pulse  shapes.  In  addition,  there  is  an  inherent  ambiguity  in 
calculating  the  phase  difference:  do  the  gamma-rays  lag  the  radio  by  0,  or  lead  by  1  —  0?  We 
ignore  these  difficulties  because  the  limited  data  do  not  warrant  a  detailed  study  here  (see 
Ulmer  1993  for  further  discussion).  We  will  remark  that  the  trend  of  increasing  offset  with 
increasing  period  predicts  that  fast  (periods  ~  50  ms)  pulsars  which  have  X-ray  emission 
dominated  by  non-thermal  processes  such  as  PSR  B0540— 69  (Seward,  Harnden,  &  Helfand 
1984;  Finley  et  al.  1993)  will  have  a  small  (^  0.1)  radio/gamma-ray  phase  offset.  It  will  be 
interesting  therefore  to  make  simultaneous  X-ray/optical  and  radio  observations  of  pulsars 
such  as  PSR  B0540— 69. 

Because  the  spectra  of  PSR  B1509— 58  from  ~  0.5  to  3.5  keV  (Seward  &  Harnden  1982; 
Seward  et  al.  1985)  and  from  ~  20  to  200  keV  (Wilson  et  al.  1993;  Matz  et  al.  1993)  appear 
to  be  well-described  by  power  laws  with  similar  spectral  indices,  the  X-ray  and  gamma-ray 
emission  is  probably  non-thermal  in  origin.  Based  on  the  apparent  energy-independence 
of  the  peak  phase  of  the  X-ray/hard  gamma-ray  light  curves  of  PSR  B1509— 58  and  of 
the  Crab  pulsar  (which  also  has  a  power-law  spectrum),  we  suggest  that  a  signature  of 
non-thermal  radiation  processes  is  the  phase  locking  of  the  ^  2  keV  emission  over  a  broad 
(^  400  keV)  energy  range.  Conversely,  the  signature  of  thermal  processes  is  the  absence 


8 


of  such  phase  locking  from  the  soft  to  the  hard  X-rays.  The  Vela  pulsar  is  an  example 
of  a  source  of  thermal  X-ray  emission  (Ogelman  et  al.  1993)  and  non-thermal  gamma-ray 
emission  (Strickman  et  al.  1993;  Kanbach  et  al.  1980)  whose  pulse  phase  varies  over  this 
energy  range. 

We  thank  the  CGRO  flight  operations  team  lead  by  B.  Breshears  for  providing  us  with  a 
stable  and  accurate  clock.  We  are  also  grateful  to  N.  Kawai  for  providing  us  with  his  data 
and  comments  on  the  Ginga  timing,  and  to  J.  Arons  and  A.  Harding  for  useful  discussions. 
We  thank  D.  Thompson  and  J.  Fierro  for  making  their  results  on  PSR  B1055— 52  available 
to  us  prior  to  publication.  This  work  was  supported  in  part  by  NASA  grant  DPR  S-10987C. 
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Table  1:  RADIO  EPHEMERIS  FOR  PSR  B1509-58 


V 

6.6371895831124  s”1 

V 

-6.76844  x  10“n  s“2 

V 

1.97  x  10-21s-3 

to  (TDB) 

48420.0 

to  geo(MJD) 

48420.000000121 

R.A.  (J2000) 

15h  13m  55f667 

Dec  (J2000) 

-59°  08'  9"42 
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FIGURE  CAPTIONS 


Figure  1.  Three  average  light  curves.  The  smooth  curve  on  the  lower  two  plots  is  the 
average  light  curve  from  K91,  normalized  and  shifted  to  fit  the  the  CGRO  data. 

Figure  2.  The  best  fit  phases  of  the  CGRO  pulses  as  a  function  of  energy.  The  K91  light 
curve  was  used  as  the  template  for  determining  the  phase.  The  triangle  is  the  data  point 
from  K91,  the  asterisks  are  the  OSSE  data,  and  the  diamonds  are  the  BATSE  data. 
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